Quasiparticle undressing in a dynamic Hubbard model: exact diagonalization study 
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Dynamic Hubbard models have been proposed as extensions of the conventional Hubbard model 
to describe the orbital relaxation that occurs upon double occupancy of an atomic orbital. These 
models give rise to pairing of holes and superconductivity in certain parameter ranges. Here we 
explore the changes in carrier effective mass and quasiparticle weight and in one- and two-particle 
spectral functions that occur in a dynamic Hubbard model upon pairing, by exact diagonalization 
of small systems. It is found that pairing is associated with lowering of effective mass and in- 
crease of quasiparticle weight, manifested in transfer of spectral weight from high to low frequencies 
in one- and two-particle spectral functions. This 'undressing' phenomenology resembles observa- 
tions in transport, photoemission and optical experiments in high Tc cuprates. This behavior is 
contrasted with that of a conventional electron-hole symmetric Holstein-like model with attractive 
on-site interaction, where pairing is associated with 'dressing' instead of 'undressing'. 



I. INTRODUCTION 

The conventional understanding of superconductivity 
starts from a normal state composed of weakly interact- 
ing 'dressed' Landau quasiparticles When the system 
goes superconducting these quasiparticles become corre- 
lated in Cooper pairs. As a function of increasing cou- 
pling strength, the Cooper pairs undergo a crossover from 
weakly bound with a long coherence length to strongly 
bound with a short coherence length. Because of the in- 
creased correlation of the quasiparticles in the Cooper 
pair their effective mass increases and their quasiparticle 
weight decreases compared to the normal state, the more 
so the stronger the coupling. In other words, the quasi- 
particle 'dressing' is larger in the superconducting than 
in the normal state. 

Instead, the theory of hole superconductivity ||-^] 
proposes a new paradigm to describe superconductiv- 
ity which is exactly opposite to what is described above. 
It also starts from a Fermi liquid of weakly interacting 
dressed quasiparticles, and Cooper pairs are also formed 
as the system goes superconducting. However, here the 
quasiparticles 'undress' and resemble more free particles 
when they are bound in the Cooper pair than when they 
are unbound in the normal state, with the effect being 
largest precisely in the strong coupling short coherence 
length regime. 

This paradoxical scenario is described by a new class 
of model Hamiltonians recently introduced to describe 
correlated electrons, 'dynamic Hubbard models'. |H,||j6) 
It is argued that these models capture an essential aspect 
of the physics of correlated electrons in atoms, molecules 
and solids, which is left out in other models like the con- 
ventional Hubbard model. Whether these new models 
describe superconductivity in any or all real materials 
remains to be established. The fundamental feature dis- 



tinguishing these from conventional models is that they 
describe 'undressing' instead of 'dressing' when carriers 
pair. The undressing should be most apparent for su- 
perconductors with high critical temperature. In sup- 
port of these models over conventional models we re- 
mark that many aspects of the phenomenology of high 
Tc cuprate superconductors, to be reviewed later, indi- 
cate that quasiparticles undress when they pair [[7|-[l2[. 
Of course there could be other unconventional models 
describing similar physics. 

In a recent paper we have begun a numerical study 
of a particular realization of a dynamic Hubbard model 
with auxiliary spin degrees of freedom. We studied 
the pair binding energy and lowering of kinetic energy 
that occurs when carriers pair as function of the param- 
eters in the model by exact diagonalization of small clus- 
ters and quantum Monte Carlo simulations, and obtained 
the approximate phase diagram of the model in one di- 
mension. The purpose of this paper is to learn more 
about the properties of this model by studying frequency- 
dependent spectral functions. The reason for doing so is 
that the properties of this simple model are likely to be 
generic for the entire class of dynamic Hubbard models 
and, as argued elsewhere representative of the prop- 
erties of correlated electrons in solids. We calculate the 
frequency dependent conductivity and the single particle 
spectral function in this model by exact diagonalization 
of small clusters, as well as the quasiparticle weight and 
effective mass, and examine their behavior as function of 
parameters. The results support the general scenario of 
'undressing' in these models indicated by more approx- 
imate treatments and qualitative arguments. Wc con- 
trast the behavior of this model with that of a conven- 
tional (electron-hole symmetric) model where 'dressing' 
rather than 'undressing' occurs upon pairing. Finally, 
we discuss the connection of our results to experimen- 
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tal observations in spectroscopic experiments in high 
materials. 



II. MODELS 

A. Dynamic Hubbard model 

The dynamic Hubbard model of interest here is defined 
by the Hamiltonian |^,|f^ 

H = Y,m - i ^[cLc+i,. + h.c] (f) 

i i.a 

where the local Hamiltonian Hi in electron representa- 
tion is 

Hi = Looa]. + gujoal + [U ~ 2guiQ(jl]ni^nii (2a) 
and in hole representation 

Hi = ujQal + gwo[2(ni| + - ^(^l 

+ [U - 2guJocrl,]ni^nii (2b) 

and al. , cr* are Pauli matrices associated with an auxiliary 
spin-f/2 degree of freedom at each site. In what follows 
the Hamiltonian in hole representation will be used. U 
is the effective on-site interaction. Briefly, the auxiliary 
spin is introduced to allow for the fact that two electrons 
on a site can be in more than one state: depending on 
the orientation of the auxiliary spin, they will experi- 
ence more or less Coulomb repulsion and in turn pay less 
or more in kinetic and single-electron potential energy, 
described by the energy of the auxiliary spin. A more 
detailed justification of the site Hamiltonian Eq. (2) to 
describe the physics of real atoms is discussed in ref . . 

In the antiadiabatic limit — * oo the effective low 
energy Hamiltonian for low hole concentration is 

-^e// = - y^[^2 + M{ni-a + n^+i,_er)](4^Ci+ier + h.C.) 
i.(T 

+U^ni^nii (3a) 



B. Electron- hole symmetric Holstein-like model 

Conventional electron-boson models involve coupling 
of a boson (e.g. phonon) degree of freedom to the elec- 
tronic charge density, and are electron-hole symmetric. 
We will contrast the behavior of the dynamic Hubbard 
model with that of a model with site Hamiltonian 

Hi = ujaal + gujocrl [n^i -I- tiii - 1] + J/o^iT^a (5) 

as a generic 'conventional' model. Some properties of this 
Hamiltonian were discussed in ref. p^ ] . The effective on- 
site interaction is 

U = Uo-2iUoiVl+9^ -I)- (6) 

The low energy effective Hamiltonian in the antiadiabatic 
limit with luq — > oo, g fixed and U fixed (i.e. Uq oo 
also) is 

Heff = -teff ^(cLq+i^ + h.c.) + UY^ Ui^n.i (7a) 

i,cr i 

with 

teff = tS' (7b) 



5 =<0|1>=<I|2 >= -(!+ ) (7c) 

V 2 yi+F 

and it gives rise to pairing only if [/ < (attractive Hub- 
bard model). For finite loq we find that the condition for 
pairing is even more stringent than U < 0. In contrast 
to the dynamic Hubbard model, here S does not become 
small as g —^ oo. 

III. SPECTRAL FUNCTIONS 
A. Optical conductivity 



t2 = SH 



At ^ tS{l - S) 
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(3b) 
(3c) 



(3d) 



which describes ground-state to ground-state transitions 
of the spin degree of freedom at each site when the holes 
hop. This Hamiltonian gives rise to pairing of two holes 
in a full band if the condition 



U g' 



At - 1+g^ 



(4) 



is satisfied. As shown in Ref [|T^ , the condition for pairing 
for finite is considerably less stringent than Eq. (4). 



We will compute the optical conductivity in these mod- 
els at zero temperature, given by 

, N I < 01 Jim > P , ^^ . ^ 

a, (c.) = ^ ^ ' F ' ^(^ - (Em - Eo)) (8) 

m 

with the current operator given by 



(9) 



It is easily shown following the steps in Maldague's 
derivation [Q4| that the sum rule 



/•oo ^ 

/ dwCTi(w) = - < 0| -T|0 > 
Jo 2 



(10) 
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holds for these models, with T the bare kinetic energy 
operator 

T^~tY,[clc,+^,a + h.c]. (11) 

When the frequency wo is not very small there is a 
natural separation of energy scales in the Hamiltonians. 
A low-lying submanifold of states in the Hilbert space 
corresponds to states where each site spin is at its site 
ground state, which is different depending on the elec- 
tronic occupation ||l^. The effective Hamiltonian Eq. 
(3) or Eq. (7) describes the coupHng between those low- 
lying states, where the electrons hop from site to site and 
the spins make ground-state to ground-state (diagonal) 
transitions. This part of the Hilbert space corresponds to 
the quasiparticle band, and the optical absorption involv- 
ing transitions between those states is the low frequency 
'intraband' part of <Ji, corresponding to the Drude part 
of the optical conductivity. We can then decompose the 
integral of the optical conductivity as 

/ dujai{uj) = / duj(Ti{uj) + / dLjai{uj) = Ai + Ah 

Jo Jo J(^m 

(12) 

where Ai is the intra-band absorption and ujm is a fre- 
quency cutoff that restricts optical transitions to the sub- 
set of intraband states. On the other hand, an optical 
sum rule also holds for the effective Hamiltonians 

/•oo 

dcual"{uj) = I < 0| - Te//|0 > (13) 

where 

Teff = - ^[^2 + At{nt-a + n.i+i-a)]{cl^Ci+la + h.C.)] 

(14) 

for the site Hamiltonian Eq. (2) and 

Teff = -teff + h.C.) (15) 

for the site Hamiltonian Eq. (5), and where crj^'^ is com- 
puted from Eq. (8) with the eigenstates and eigenval- 
ues of the effective Hamiltonians. Because the low-lying 
spectrum and eigenstates of the full Hamiltonians coin- 
cide with those of the effective Hamiltonians we can write 
Eq. (13) as a 'partial' conductivity sum rule for the low- 
lying eigenstates and eigenvalues of the full Hamiltonian 

/ di^ai{uj)^ Ai^ - <0\-Teff\Q> (16) 
Jo ^ 

The ground state |0 > on the right side of Eq. (16) is the 
ground state of the effective Hamiltonian Eq. (3) or Eq. 
(7). 



For the effective Hamiltonian of the dynamic Hubbard 
model Eq. (3) it can be seen within BCS theory [|l5[ 
as well as from the exact solution in the dilute limit | [l^ ] 
that when pairing occurs the right-hand side of Eq. (16) 
increases. This extra spectral weight signals a lowering of 
kinetic energy and effective mass reduction in the quasi- 
particle band. In a real physical system the total inte- 
grated optical spectral weight is conserved, so that any 
extra spectral weight at low frequency has to come at the 
expense of spectral weight in another frequency range. In 
tight binding models however the total optical spectral 
weight is not conserved, because the current and kinetic 
energy operators do not describe transitions to states in 
other bands. As a consequence, the total integral Eq. 
(12) does not remain constant but also increases upon 
pairing in our model. Nevertheless, we will see in the 
numerical results that the dominant effect of pairing in 
the dynamic Hubbard model is a large increase in the low 
frequency 'intra-band' special weight Ai and an overall 
shift in optical spectral weight from higher to lower fre- 
quencies. 

In contrast, we will see that in the electron- hole sym- 
metric model, the behavior is exactly opposite: in the 
paired state the effective mass of the carriers increases 
hence optical spectral weight is transfered from the quasi- 
particle band to higher frequencies. 

B. One-particle spectral function 

The changes in the optical conductivity upon pairing 
are intimately related to changes in the one-particle spec- 
tral functions. We consider here the spectral function for 
hole destruction in a system of ?i + 1 holes, defined as 

A„+i,„(c^) - ^ I < Z„|cfc.|0„+i > f8{Lo - [El' - E^)) 
I 

(17a) 

as well as the spectral function for hole creation in a 
system of n holes: 

A„,„+i(c.) - 5] I < ;n+i|cL|0„ > I'Siu - {E-+^ - E-)) 
I 

(17b) 

Here, |Z„ > denotes the l—th excited state of the system 
with n holes. In the exact diagonalization calculation we 
will study clusters with one hole and two holes of opposite 
spin. The relevant momentum is then fc = 0. 

Consider these functions for a single site in the dy- 
namic Hubbard model. They are given by 

A2iiu) ^ Auiu) = S{u) (18a) 

Awiiu) = Aoi(cj) = S^6{iu) + (1 - S^)Siu; - 2uoV'^ + 9^) 

(18b) 
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with S given by Eq. (3d). The lo = part corresponds 
to the quasiparticle contribution, and its coefficient is the 
quasiparticle weight z. It can be seen that for the single 
site, z = 1 for the spectral function involving one and 
two holes and z < 1 for the spectral function involving 
zero and one hole. In an extended system, when z < 1 
there will be an 'incoherent' contribution to the spectral 
function represented in the site by the second term in Eq. 
(18b). Hence we expect in an extended system that in 
the unpaired state the quasiparticle weight will be small, 
and that when pairing occurs spectral weight in the sin- 
gle particle spectral function will be transfered from the 
high frequency incoherent part to the a; = quasiparti- 
cle peak. Correspondingly, optical transitions involving 
transitions between singly and doubly occupied sites will 
have larger low energy spectral weight than those involv- 
ing empty and singly occupied sites, and optical spectral 
weight should be transfered from high frequencies to low 
frequencies when pairing occurs. 

In contrast, in the electron- hole symmetric model the 
spectral functions for creation of a hole in a singly occu- 
pied site and destruction of a hole in a singly occupied 
site are both given by 

Ai2(w) = ^io(c^) - S^S{lu) + (1 - S^)6iiJ - 2coov/l + .g2) 

(19a) 

and the spectral functions for destruction of a hole in a 
doubly occupied site and creation of a hole in an empty 
site are given by 

A2i{uj) = Aoi(w) = S^5{uj) + (1 - S^)S{uj - 2wo) (19b) 

where S is given by Eq. (7c). It can be seen that the 
quasiparticle weight is the same for all these spectral 
functions (z = S*^), whether the site is initially empty, 
singly or doubly occupied. Hence we cannot extract any 
conclusions about changes in the quasiparticle weight 
upon pairing from 'single site' physics in this model. In 
an extended system however because of the more corre- 
lated nature of the wavefunction in the paired state we 
will see that pairing is associated with decrease of the 
quasiparticle weight, i.e. increased 'dressing', in contrast 
to the behavior in the dynamic Hubbard model. Simi- 
larly one would expect transfer of optical spectral weight 
from low frequencies to high frequencies upon pairing in 
this model, in contrast to the behavior in the dynamic 
Hubbard model. 



IV. NUMERICAL RESULTS 

We diagonalize exactly the Hamiltonian Eq. (1) with 
site Hamiltonians Eq. (2b) (dynamic Hubbard model) 
and Eq. (6) (electron-hole symmetric model) for four 
sites in the subspaces with zero-, one- and two- hole oc- 
cupation. We choose units so that t — 1. For a finite 
chain, the optical sum rule Eq. (10) holds if free ends 



boundary conditions are use, but not if periodic bound- 
ary conditions are used; in the latter case, an extra term 
proportional to a zero-frequency (5-function is needed to 
satisfy the sum rule (l^. We will use free ends bound- 
ary conditions to calculate the optical conductivity, so 
that the total optical spectral weight is obtained in the 
numerical calculation. The lowest frequency peak then 
occurs at a finite frequency ujd ('Drude 'precursor') |l^] 
that goes to zero as the size of the system increases. 

A. Dynamic Hubbard model 

To obtain a clear separation of energy scales that illus- 
trates clearly the physics of the model we choose a strong 
coupling case, with g = 3. The single hole effective hop- 
ping Eq. (3b) in the antiadiabatic limit is then t2 — 0.1, 
and At — 0.216. The site energies for n holes (excluding 
the U term) are given by 

e{n) = -e{n) = -uJoVlT^ (20) 

with e{n),e{n) the energies of the site ground state and 
excited state. Hence the difference e(n) — e(n) = 6.3a;o 
is much larger than the effective bandwidth for the holes 
except for very small luq. Even so, for the case of 2 holes 
the effective Hamiltonian Eq. (3) is not an accurate rep- 
resentation for finite frequencies because of the contribu- 
tions from 'vertical transitions' as discussed in ref. p3| ; 
the effective interaction 

Ueff ^ E{2)+E{0)~2E{1) (21) 

{E{n) =energy with n holes) is considerably more attrac- 
tive for finite frequency than in the antiadiabatic limit 
ciJo ^ oo where the effective Hamiltonian Eq. (3) is valid. 
We will use loq — 1 here. 

Figure 1 shows the effective interaction versus the on- 
site repulsion U . For luq = 1 \t is attractive for U < 5.9, 
while in the antiadiabatic limit it is attractive only for 
U < 3.6. 

It is interesting to consider the structure of energy lev- 
els, shown in Figure 2. We show the energy levels for a 
single hole in the system and for two holes in the cases 
U = 8 and U — 0. The effective interaction for these two 
cases is Ue/f — 0.056 and Ueff = —0.715 respectively. 
There is a clear separation between low-lying energy lev- 
els, described by an effective 'intraband' Hamiltonian and 
higher-lying levels where the spin degrees of freedom are 
in excited states. In the single particle case there are 
4 low-lying states, in the two-particle case there are 16 
low-lying states for U ~ , ior U — 8 there are only 
12 because 4 'intraband' states are pushed high in en- 
ergy due to the strong on-site repulsion. Note also that 
the energy range of the low-lying states is considerably 
larger for the case U = than for = 8. This reflects 
the bandwidth expansion that occurs upon pairing in this 
model. 



4 



Figure 3 shows the frequency-dependent conductivity 
for on-site repulsion U ~ 8 (solid line) and for U — 
(dashed line). For U = 8 the holes are not bound since 
Ueff > 0, while for U = the holes are bound with 
Ueff = —0.715. The low frequency 'intraband' conduc- 
tivity is seen to increase substantially when the holes are 
paired. Furthermore, the optical absorption at interme- 
diate frequencies (wq ~ 6) increases upon pairing, while 
the optical absorption at the highest frequencies {lo > 10) 
decreases upon pairing. This shows that optical spectral 
weight is transfered from high to low frequencies when 
pairing occurs. Interestingly, not only the 'intraband' 
optical spectral weight increases but also spectral weight 
at intermediate frequencies. This will be discussed fur- 
ther in the next section. Note that changes in optical 
spectral weight occur at energies much higher than the 
scale of the pairing energy. |C/e//| ^ 0.7. 

In Figure 4a we show the optical spectral weight asso- 
ciated with the intra-band optical transitions Ai as well 
as the total optical spectral weight Ai + Ah versus on- 
site repulsion U. As mentioned earlier the model does 
not conserve total optical spectral weight, hence Ai + A^ 
is not a constant; nevertheless, the intraband spectral 
weight increases faster than the total spectral weight as 
pairing occurs, and the ratio Ai/A^ increases as U de- 
creases as shown in Figure 4b. Hence the model does 
describe a transfer of optical spectral weight from high 
frequencies to low frequencies as pairing occnirs. 

Figure 4a also shows the optical spectral weight in the 
antiadiabatic limit u>o ^ oo. As increases the high 
frequency optical spectral weight moves to higher fre- 
quencies and decreases in amplitude, and vanishes for 
Wo = 00. However the low frequency cri(ti;) shows very 
little dependence on ojq and its total weight is almost the 
same for loq = 1 and wq = oo as seen in Figure 4a. 

The 'undressing' process that occurs upon pairing is 
also clearly seen in the one-particle spectral functions. 
Figure 5 shows the spectral functions for single hole de- 
struction in the system with one hole, Aiq^lv), and in 
the system with two holes, ^21(1^)- The lowest frequency 
peaks are actually ^-functions at zero energy, their weight 
is the quasiparticle weight. For the system with a sin- 
gle hole, the quasiparticle weight is z = 0.174. This is 
larger than the quasiparticle weight for a single site for 
these parameters (z = 0.1); in the extended system, fi- 
nite ajQ gives rise to a larger quasiparticle weight due to 
retardation. The quasiparticle weight for the two-hole 
system and large on-site U is very similar to the single 
hole spectral weight, z = 0.176. When U decreases and 
a pair is formed, the quasiparticle weight increases, to 
z = 0.422 when [7 = 0. As seen in Figure 5b, spectral 
weight is transfered from the incoherent region of the 
spectral function at energies around a; 6 to the quasi- 
particle peak. This energy range corresponds to states 
where one of the background spins in the system is in 
a final excited state. There is also a very small hut fi- 
nite spectral weight in the energy region around w ~ 12, 
where two of the background spins are in excited states 



in the final state. 

Figure 6 shows the spectral function for hole creation 
in a system of zero and one hole, Aqi^lo) and Ai2(w). 
Unlike the spectral functions for hole destruction they 
satisfy the sum rule 

POO 

/ cLjAn,n+l{u^) = 1 (22) 

Jo 

(n = or 1) because they have no negative frequency 
component. The quasiparticle weights extracted from 
these spectral functions are of course the same as those 
obtained from Figure 5. It is interesting to note that in 
the intra-band energy range {lo < 2) there is now an in- 
coherent contribution to the spectral function A12 (w) for 
parameters where the effective interaction is attractive. 
The 'intraband' part of the ground state pair wavefunc- 
tion is of the form 

|*>=^/(fc)4/_,j0> (23) 

k 

and it has finite overlap with states other than cj^cjj^ 10 >. 
Still the weight of the coherent part of the spectral func- 
tion, i.e. the zero frequency peak (quasiparticle weight) 
increases strongly as the pair is formed. 

The behavior of quasiparticle weights as function of on- 
site repulsion is shown in Figure 7. For large repulsive 
U, the quasiparticle weight for a hole is the same for the 
system with two holes (dashed line) and with one hole 
(full line); as the on-site repulsion decreases and the pair 
forms, the quasiparticle weight in the system with two 
holes increasing, indicating that 'undressing' occurs. We 
also show the corresponding results in the w = 00 limit, 
which show similar behavior; the quasiparticle weights 
in this case are smaller than for finite frequency. For the 
single hole, the quasiparticle weigth in the antiadiabatic 
limit is the same as for the single site, z = 5^. 

In summary, these results show that in the dynamic 
Hubbard model there is 'undressing' when pairing oc- 
curs: spectral weight in both one- and two-particle spec- 
tral functions is transfered from high to low frequencies, 
the quasiparticle weight increases and the kinetic energy 
decreases when pairing occurs. We next discuss the sit- 
uation for the electron-hole symmetric model. 

B. Electron- hole symmetric Holstein-like model 

We consider the electron-hole symmetric (e-h sym) 
model Eq. (5), for parameters g = 3 and wq = 1- The 
site energies for n holes in this case are 

e(n) = -e(n) = -tooVTTg^ (24a) 
for n = and n = 2, and 

e(n) = -e(n) = -ujq (24b) 
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for n = 1. Because the excitation energy of the back- 
ground spin is lower for the singly occupied site here 
there is not such a clear separation of energy scales for 
this model as for the dynamic Hubbard model for the 
parameters used. 

Figure 8 shows the effective interaction between two 
holes in the four-site cluster versus on-site repulsion, as 
well as the results in the antiadiabatic limit loq — oo. 
Note that retardation is detrimental to pairing in this 
model. For wq = 1, the effective interaction is attractive 
for on-site repulsion U smaller than —0.4. 

The optical conductivity for two holes in the cluster is 
shown in Figure 9, for on-site repulsion U = 2, where the 
effective interaction is repulsive, and for [/ — —2 where 
it is attractive. In this case a transfer of optical spectral 
weight from low to high frequencies occurs upon pairing, 
opposite to the behavior in the dynamic Hubbard model 
seen in Figure 3. Note also that the separation between 
low and high energy regions is less clear here than in the 
previous case, as expected. 

The dependence of optical spectral weights on the on- 
site repulsion is shown in Figure 10a. The intra-band 
spectral weight decreases by a factor of 2 as pairs form; 
this indicates that the carriers are more heavily dressed 
and have larger effective mass in the paired state, as one 
would expect; in the lo oo limit this model becomes 
equivalent to an attractive Hubbard model, where the 
pair mobility is always smaller than the single particle 
mobility p8[ ; in particular in the strong coupling limit 
the pair hopping is tp — 2t^ /\U\ much smaller than the 
single particle hopping for large on-site attraction. This 
model then describes a transfer of optical spectral weight 
from low to high frequencies when pairing occurs; the ra- 
tio of intra-band to inter-band optical spectral weights 
versus on-site U is shown in Fig. 10b, qualitatively dif- 
ferent from the corresponding results for the dynamic 
Hubbard model Figure 4b. 

The single particle spectral functions for hole destruc- 
tion are shown in Figure 1 1 . The quasiparticle weight for 
the system with a single hole is z = 0.89, for the system 
with two holes it is z = 0.79 for the case U — 2 with 
unpaired holes, and it decreases to z — 0.64 for U = —2 
when the holes are paired. Note that this is not a single- 
site effect, as for a single site the spectral weight for one 
and two holes is identical due to electron-hole symmetry 
(Eq. (19)). Figure lib shows that spectral weight from 
the quasiparticle peak is transfered to higher frequencies 
around w ^ 3 when pairs form. Similarly, Figure 12 
shows the single particle spectral function for hole cre- 
ation in the system with zero and with one hole. For 
the creation of a hole in the empty cluster the incoher- 
ent part is here at lower frequencies due to the lower 
excitation energy of the singly-occupied site, Eq. (24b). 
Again, figure 12b shows that as pairing occurs spectral 
weight is transfered from the quasiparticle peak to higher 
frequencies. 

Figure 13 shows the variation of quasiparticle weight 
versus on-site repulsion in this model. Because for the 



single site the quasiparticle weight is independent of hole 
occupation the variation here is less than in the dynamic 
Hubbard model (Fig. 7). The quasiparticle weight in the 
system with 2 holes decreases as the pair formed, indi- 
cating that the quasiparticle is more heavily dressed in 
the paired compared to the unpaired state, qualitatively 
different to the situation in the dynamic Hubbard model. 

Figure 13 also shows the quasiparticle weights in this 
model in the antiadiabatic limit, where it is equivalent to 
the Hubbard model. Note that the quasiparticle weights 
for the system with 1 and 2 holes coincide for U = 0, 
where the Hamiltonian describes non- interacting holes, 
and are both given by the site value z = S'^ = 0.658; both 
for repulsive and attractive U the quasiparticle weight is 
smaller in the two-hole system due to 'intra-band' hole- 
hole interaction. The suppression of quasiparticle weight 
however is larger for negative than for positive U . 



C. Summary 

The qualitatively different behavior of the two models 
considered is summarized in Figure 14. We plot the effec- 
tive mass enhancement and the quasiparticle weight as 
function of the effective interaction Ueff in both models. 
The ratio of effective masses for the single particle and 
the pair, m*/m*, is obtained from the ratio of the intra- 
band optical spectral spectral weight for the system with 
two holes and twice the intra-band optical spectral weight 
for the system with one hole. Figure 14a shows that in 
the dynamic Hubbard model the pair becomes increas- 
ingly lighter than the single particle as the interaction 
becomes more attractive, while in the electron- hole sym- 
metric model the pair becomes increasingly heavier as 
the interaction becomes more attractive. Similarly, the 
quasiparticle weight increases in the dynamic Hubbard 
model as the pair is formed, while it decreases in the 
electron-hole symmetric model. In summary, the system 
becomes more coherent in the paired state in the dynamic 
Hubbard model and more incoherent in the electron-hole 
symmetric model. 

It is interesting to note that the quasiparticle weight 
and the effective mass change by approximately the same 
factor in the dynamic Hubbard model in the parameter 
range considered, 2.5 and 2.4 respectively. This is what 
one would expect if the physics of the model is dominated 
by 'single site' physics, where the quasiparticle weight 
and effective hopping and hence effective mass are both 
determined by the single site overlap matrix element S. 
More generally, in a many body system the exact single 
particle Green's function is given by p9[ 



Zk 



w - ek - T.{k,uj) uj - Ek 



-+G'{k,uj) (25) 



with E(fc,w) the self-energy, ik the quasiparticle energy 
and G" the incoherent part of the Green's function. The 
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quasiparticle weight Zk and effective mass enhancement 
are given by 



d 

Zk = (1 - — Sre(fc,Cj))"^ 
OLU 



m dik 
m* oek 



) 



(26a) 



(26b) 



(Ere = real part of E) so that if the self-energy is mo- 
mentum independent the quasiparticle weight and effec- 
tive mass renormalization coincide. Hence our results 
indicate that this is approximately the situation in the 
dynamic Hubbard model, and suggest that dynamical 
mean field theory | pO| , which assumes a momentum- 
independent self-energy, should be a useful approach to 
study these models. In contrast, the results for the 
electron-hole symmetric model in the parameter range 
considered yield an effective mass changing by a factor 
1.9 with the quasiparticle weight changing by a factor of 
only 1.2 , suggesting that for this model the momentum- 
dependence of the self-energy is substantial. 



D. Finite temperatures 

We have also studied the behavior of the optical con- 
ductivity in the dynamic Hubbard model at finite tem- 
peratures, given by 



En 



En 



-|<0|J|m>|2x (27) 

5{lO - [Em ~ En)) 



with Z = g-/3s„ ^Yie partition function. Figure 15 
shows 0-1(0;) for a case where the effective interaction is 
attractive, C4// = —0.715. The intra-band part of the 
optical absorption is rapidly suppressed as the tempera- 
ture increases, and the peak at intermediate frequencies 
is also suppressed. Unfortunately because tight bind- 
ing models do not satisfy the optical sum rule the to- 
tal optical spectral weight in the model is not conserved 
as the temperature changes. Nevertheless it is interest- 
ing to note that as T increases optical spectral weight 
is transfered to the very high frequency region oj > 10. 
This change in high frequency spectral weight occurs on 
a temperature scale (T ~ 1) that is related to the energy 
scale of the pairing energy (C/g//) and unrelated to the 
energy scale of the frequency where the optical absorp- 
tion change occurs. 



V. RELATION WITH EXPERIMENTS 

The dynamic Hubbard model describes a coupling of 
electrons to a background degree of freedom (the auxil- 
iary spin) that only exists when the site is doubly occu- 



pied by electrons. Hence its effect decreases with increas- 
ing local hole concentration, on the average the coupling 
constant is 



A(n) = (2 - n) 



(28) 



with n the average hole concentration per site. This leads 
to a phenomenology whereby hole carriers 'undress' in 
the presence of other hole carriers. The undressing man- 
ifests itself in transfer of optical spectral weight from high 
to low frequencies and in transfer of one-particle spec- 
tral weight from the incoherent (high frequency) to the 
coherent (low frequency) region; these spectral weight 
transfers lead to increase in the quasiparticle weight and 
decrease in the quasiparticle effective mass. Because a 
hole comes close to another hole both when the hole con- 
centration is increased by doping in the normal state as 
well as when carriers pair, undressing will occur both for 
increasing carrier concentration and for decreasing tem- 
perature. 

There is substantial evidence for such phenomenology 
in the high Tc cuprates. Johnson et al |9| extract from 
photoemission experiments in the normal state a coupling 
constant that decreases continuously as the hole doping 
increases, as well as a mass enhancement that decreases 
with hole doping. In the overdoped regime, Yusof et al 
pl| find evidence from photoemission for the existence 
of quasiparticles in the normal state, which appear to be 
absent in the underdoped regime. This is consistent with 
the phenomenology of the dynamic Hubbard model in a 
strong coupling regime where the quasiparticle weight for 
low hole density would be small enough to be unobserv- 
able. The quasiparticle weight in the dynamic Hubbard 
model as function of hole concentration is approximately 
given by 



zin)^S\l + ^T)^ 



S 



- 1 



(29a) 



(29b) 



If the coupling strength g is large S will be small, z be- 
comes very small in the underdoped regime, and at the 
same time the 'undressing parameter' T becomes large, 
hence quasiparticles undress rapidly with increasing n. 

Ando et al |^ extract from transport measurements 
a hole mobility that increases monotonically with hole 
doping, consistent with the behavior predicted by the 
dynamic Hubbard model and its low energy effective 
Hamiltonian Eq. (3), which leads to an effective density- 
dependent hole hopping 



ieff{n) = t2 + nAt 



(30) 



which is equivalent for low n to the more fundamental 
relation 



teffin) = tz{n) 



(31) 
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describing the fact that the quasiparticle effective mass 
is inversely proportional to the quasiparticle weight as 
expressed by Eq. (26b) when the self-energy has no mo- 
mentum dependence. 

Ding et al Q make a compelling description of the 
undressing phenomenology of high Tc cuprates: from 
their photoemission data they extract a quasiparticle 
weight that emerges from an incoherent background as 
the temperature is lowered and the system becomes su- 
perconducting. This is of course consistent with the phe- 
nomenology of the dynamic Hubbard model in the strong 
coupling regime where z will increase strongly when pair- 
ing occurs, as seen in Figs. 5 and 7. Furthermore Ding 
et al find that z increases as the doping increases, again 
consistent with the behavior expected in the dynamic 
Hubbard model. 

In the finite cluster calculations reported here for the 
optical conductivity we cannot distinguish whether the 
optical spectral weight transfered to low frequencies when 
pairing occurs goes into the zero-frequency (5— function or 
to finite frequencies. However the effective Hamiltonian 
for the dynamic Hubbard model clearly describes transfer 
of spectral weight to the (5— function as the pairing ampli- 
tude develops ||T^,|I^, since the average effective kinetic 
energy 

^e// = - Y^'^ff^'^^ < (^\aCi+la + h.C.) > (32) 

-2At{< c\^c\^ >< Ci„Ci+i„ > +h.c.] 

has a contribution from anomalous expectation values, 
while the intra-band optical absorption is unchanged ex- 
cept for the depletion due to the opening of the supercon- 
ducting energy gap |2^. This indicates that the optical 
spectral weight transfered into the intra-band region goes 
into the zero-frequency (5— function, and will lead to an 
apparent violation of the Ferrell-Glover-Tinkham optical 
sum rule |ll5|| . This violation has been observed experi- 
mentally, both for in-plane as well as c-axis light polar- 
ization by Santander-Cyro et al and by Basov et al 
10 respectively. The decrease in high frequency spectral 
weight predicted by the dynamic Hubbard model (range 
above w ~ 10 in Figs. 3 and 15) is also consistent with 
recent experimental observations by Molegraaf et al , 
who report a decrease in optical spectral weight in the 
energy range between 1.25eV and 2.5eV which is trans- 
fered to lower frequencies. This is also consistent with 
earlier observations by Fugol et al p3|l . 

Note that our results for the dynamic Hubbard model 
predict also an increase in optical spectral weight at inter- 
mediate frequencies, w ~ 6 when the interaction becomes 
attractive (Fig. 3) or when the temperature decreases 
(Fig. 15), well above the 'intra-band' frequency range. 
We expect this to be a generic feature of these models, 
describing optical transitions where a hole is transfered 
to a nearest neighbor site already occupied by another 
hole, leaving the auxiliary spin behind in an excited state. 
Such a 'vertical' transition will have a large weight in the 



paired state (proportional to 1 rather than to S if the 
neighboring site is unoccupied) and enhances the optical 
absorption at frequencies corresponding to the excitation 
energy of a single site. Given the correspondence dis- 
cussed above between the range w > 10 in our example 
and the visible frequency range, the intermediate region 

~ 6 would correspond to the mid-infrared range of fre- 
quencies in the cuprates. Indeed, Gao et al |2J] report 
observation of extra optical spectral weight appearing be- 
low Tc in the mid-infrared region. 

Table I summarizes the matrix elements for the various 
optical transitions in the dynamic Hubbard model con- 
sidered here for the case of a single hole versus the case of 
two holes on neighboring sites. The weight ratios given 
in the last column summarize the expected qualitative 
behavior of spectral weight changes when hole doping in- 
creases in the normal state or when the temperature is 
lowered and the system goes superconducting. We ex- 
pect the frequency range of the second and third rows to 
correspond to mid-infrared and that of the fourth row to 
visible frequencies in the cuprates. Because S+l/S > 2S 
for any S, the optical spectral weight in the mid-infrared 
range should always increase upon pairing. 

VI. DISCUSSION 

We have studied in this paper the behavior of spec- 
tral functions in a dynamic Hubbard model by exact 
diagonalization. The results obtained are exact for the 
small cluster studied. They support a scenario for the 
physics of this class of quantum many-body systems 
that is expected from qualitative arguments and approx- 
imate treatments. Namely, that quasiparticles 'undress' 
in these models when the local hole concentration in- 
creases, which occurs both when holes are added to the 
system (doping) and when holes pair and form Cooper 
pairs. This scenario is obtained from the exact calcula- 
tions in this paper without uncontrolled approximations. 

The physics of dynamic Hubbard models is especially 
transparent in the antiadiabatic limit, where the effective 
Hamiltonian is a Hubbard model with correlated hop- 
ping. The results of ref. ||l^ as well as the results in 
this paper indicate that the 'intra-band' physics of the 
model described by the antiadiabatic limit — > oo re- 
mains essentially unchanged with ujq decreasing to rather 
small values. In addition to the case luq = I discussed in 
this paper we also studied the model for luq = 0.5 and 
obtained qualitatively similar results. What changes as 
luq decreases is that the energy scale of non-intraband 
excitations decreases up to a point where there is no 
longer a clear separation between 'intra-band' and 'non- 
intraband' regions. Nevertheless the low energy physics 
remains essentially unchanged. 

In a more realistic description of a real system there 
will presumably be a set of excitation energies iVi describ- 
ing the bosonic excitations that dress the hole quasipar- 
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tides. Still we do not expect the physics in such a case to 
be qualitatively different. These excitations would cor- 
respond to local electronic excitations with scale up to 
eV's. Such energy scales are consistent with experimen- 
tal observations in cuprates that optical spectral weight 
in the visible range is transfered to low frequencies both 
when the system goes superconducting p2] a s well as 
when it is hole-doped in the normal state [p5| . Our re- 
sults in this paper as well as earlier results |[26| demon- 
strate that dynamic Hubbard models naturally describe 
transfer of spectral weight from high energies unrelated 
to the scale of superconductivity down to low frequen- 
cies when superconductivity sets in, as observed p2| , p3[ , 
hence they provide a natural explanation for the origin of 
the high energy scale observed in the optical experiments 
of Molegraaf et al . The detailed description of such 
apparently counterintuitive physics remains a challenge 
for other proposed descriptions of the physics of high tem- 
perature superconductors that also propose that super- 
conductivity is driven by kinetic energy lowering |^,^ . 

The dynamic Hubbard model and the electron-hole 
symmetric model discussed in this paper are represen- 
tative of two classes of model Hamiltonians, of which 
there are many different realizations. In particular, it is 
not essential that the coupling of the boson in a model 
in the class of dynamic Hubbard models be to the on- 
site double occupancy; a model with coupling only to 
the on-site charge density will also belong to this class if 
it is not electron- hole symmetric Q. Also the auxiliary 
boson may be an oscillator rather than a spin, or the 
model could have only electronic degrees of freedom [^. 
What distinguishes these two classes of models is what is 
the driving energetics for pairing: in dynamic Hubbard 
models pairing is kinetic energy driven, and the potential 
energy increases upon pairing, and the opposite is true 
in the other class of models, which may be termed 'con- 
ventional' or 'electron- hole symmetric'. The conventional 
electron-phonon models used to describe conventional su- 
perconductors belong to this second class of models. We 
believe that these two classes of models represent very 
general paradigms. In the class of dynamic Hubbard 
models the 'undressing' of quasiparticles is essential to 
lead to kinetic energy lowering; instead, in the 'conven- 
tional' class of models the dressing of the quasiparticles 
may remain unchanged upon pairing if the coupling is 
weak, or increase in a strong coupling regime. 

Dynamic Hubbard models and their low energy effec- 
tive Hamiltonians can describe superconductivity over 
the entire range of coupling strengths. The physics is 
determined by the scale of excitation energies ujo and by 
the strength of the couping g or equivalently the magni- 
tude of the 'undressing parameter' T (Eq. (28b). As we 
have seen in this paper and in ref. |13| the low energy 
physics is not strongly dependent on the scale ujq. The 
magnitude of Tc and the superconducting gap is mainly 
determined by the strength of the dimcnsionless param- 
eters g or T and the single electron hopping parameter t, 
as well of course as competing Coulomb repulsions such 



as U. These parameters cannot however be tuned sepa- 
rately at will, in a real system they are all closely interde- 
pendent and determined principally by the ionic charge 
Z as discussed in refs. For increasing g, Tc becomes 

large, the coherence length in the superconducting state 
becomes short and the system becomes incoherent in the 
normal state for low hole concentration as the magni- 
tude of S and the quasiparticle weight decrease; the 'un- 
dressing' phenomenology becomes particularly apparent 
in this regime. For small g , Tc becomes small, the co- 
herence length can become thousands of lattice spacings 
and the normal state becomes coherent; in this regime, 
even though it is still the same 'undressing physics' that 
drives the transition to superconductivity, the anoma- 
lous spectral weight transfers signaling undressing will 
become almost invisible. We suggest that it is possible 
that this same physical mechanism can describe the su- 
perconducting phenomenology of materials as distinct as 
YBa2Cu307-s, MgB2, and Al M. 
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FIG. 1. Effective interaction Eq. (21) versus on-site re- 
pulsion U for the dynamic Hubbard model with g — 3 and 
LJo = 1 (full line) and in the antiadiabatic limit uio = oo 
(dashed line). The same parameters {g = 3,uio = 1) are used 
in the following figures. 



FIG. 7. Quasiparticle weights versus U for the system with 
two holes (dashed line) and with one hole (full line). The 
corresponding quasiparticle weights in the antiadiabatic limit 
are shown as the dash-dotted and dotted lines respectively. 



FIG. 8. Effective interaction Eq. (21) versus on-site repul- 
sion U for the electron-hole symmetric model Eq. (5), with 
g — 3 and uio = 1 (full line) and in the antiadiabatic limit 
LtJo = oo (dashed line). The same parameters (g = 3,tJo = 1) 
are used in the following figures. 



FIG. 9. Frequency-dependent conductivity for the elec- 
tron-hole symmetric model with two holes and on-site repul- 
sion (7 = 2 and (7 = -2. 



FIG. 2. Energy levels for the dynamic Hubbard model. 
The energy levels for a single hole in the four-site cluster and 
for two holes with on-site repulsion U — 8 and U = are 
shown. 

FIG. 3. Frequency-dependent conductivity for dynamic 
Hubbard model with two holes and on-site repulsion U = 8 
and U — 0. The 5-functions in Eq. (8) are broadened 
to lorentzians with width P = 0.5. The lowest frequency 
5— function at frequency uid ('Drude precursor') is shifted to 
LJ — and represented by a Drude form (semi-lorentzian) with 
width r = 0.5. 



FIG. 10. (a) Kinetic energies for electron-hole symmetric 
model versus on-site U. The cutoff frequency to define the 
'intraband' spectral weight Ai is uirn ~ 2. The dash-dotted 
line gives the results for twice the total optical spectral weight 
Ai + Ah for a single hole in the cluster, and the dotted line 
gives the corresponding intra-band value, (b) Ratio of in- 
tra-band to inter-band optical spectral weight versus U. 

FIG. 11. One-particle spectral function for hole destruction 
in the electron-hole symmetric model with one hole (a) and 
with two holes (b). 



FIG. 4. (a) Kinetic energies for dynamic Hubbard model 
versus on-site U, obtained from integration of ai{uj) from Eqs. 
(12) and (16). The cutoff frequency to define the 'intraband' 
spectral weight Ai is LOm = 2. The dash-dotted line gives the 
results for twice the total optical spectral weight Ai + Ah 
for a single hole in the cluster, which coincides with the re- 
sults for two holes in the cluster for large U. The dotted 
line gives twice the intra-band optical spectral weight for one 
hole, which is approximately equal to the intra-band spectral 
weight for two holes in the cluster for large U. The intra-band 
spectral weights for two holes are shown both for = 1 and 
in the antiadiabatic limit loq = ex; for one hole in the antia- 
diabatic limit the results are indistinguishable from the value 
for LOQ — 1. (b) Ratio of intra-band to inter-band optical 
spectral weight versus U for luq = 1. 

FIG. 5. One-particle spectral function for hole destruction 
in the system with one hole (a) and with two holes (b). The 
S— functions are broadened to lorentzians with width F — 0.1. 
The spectral weight at frequencies above = 10 is very small 
and is amplified in the figure by a factor 50. 



FIG. 12. Same as figure 11 for hole creation in the system 
with zero holes (a) and with one hole (b). 



FIG. 13. Quasiparticle weights versus U for the elec- 
tron-hole symmetric model with two holes (dashed line) 
and with one hole (full line). The corresponding quasipar- 
ticle weights in the antiadiabatic limit are shown as the 
dash-dotted and dotted lines respectively. 

FIG. 14. (a) Ratio of effective mass of a single hole to the 
effective mass of a hole in a pair, calculated from the ratio of 
intra-band kinetic energies, versus effective interaction Ueff 
(Eq. (21)) for both models. For the dynamic Hubbard model 
this ratio increases as Uef / decreases and the pair forms, while 
for the electron-hole symmetric model it decreases, (b) Ra- 
tio of quasiparticle weights of a hole in the system with two 
holes to the quasiparticle weight of the single hole, versus ef- 
fective interaction. In the dynamic Hubbard model this ratio 
increases as Ueff decreases and the pair forms, in the elec- 
tron-hole symmmetric model it decreases. 



FIG. 6. Same as figure 5 for hole creation in the system 
with zero holes (a) and with one hole (b). 
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FIG. 15. Optical conductivity of the dynamic Hubbard 
model for two holes and on-site repulsion ?7 = for vari- 
ous temperatures (in units of the bare hopping t). Note the 
transfer of spectral weight from high frequencies (w > 10) to 
low frequencies as the temperature is lowered. 

TABLE I. Optical transitions, possible states of auxiliary spin. The ground state and excited state of the auxiliary spin at a 
site with n holes are denoted by \n > and \n > respectively. The four left columns correspond to transitions involving a single 
hole hopping between neighboring sites, | t> |0 >^ |0 > | t>; the four right columns correspond to transitions involving two 
holes at neighboring sites, | |> | i>— > |0 > | tJ.>. 
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1 T> |o > 


final state 
|0 > 1 T> 


weight 


energy 


initial state 
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|1>|0> 
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Figure 13 
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